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— ~ Objectives 

The development of a composite wing box section using a higher order-theory is 
proposed for accurate and efficient estimation of both static and dynamic responses. The 
theory includes the effect of through-the-thickness transverse shear deformations which is 
important in laminated composites and is ignored in the classical approach. The box beam 
analysis is integrated with an aeroelastic analysis to investigate the effect of composite 
tailoring using a formal design optimization technique. A hybrid optimization procedure is 
proposed for addressing both continuous and discrete design variables. 

Accomplishment 

In recent years, aeroelastic tailoring has received considerable attention as a means 
to improve aeroelastic performance through directional stiffening. Fibrous composite 
materials can further enhance aeroelastic tailoring capabilities by utilizing their unique 
stiffness and strength properties. Accurate and efficient prediction of the structural 
response is very important in the investigation of aeroelastic tailoring using composite 
structures. The analysis of aircraft wings can be done either through a detailed 
investigation of the wing sections comprising spars, webs, ribs etc., or through the 
representation of the load carrying member by reduced composite box beam models. The 
detailed analysis is computationally very expensive and is often impractical in design 
optimization and/or trade-off studies. Several approaches addressing composite box beam 
modeling have been proposed over the last few years [1-2]. All of these models have 
several limitations. In some of these work, classical laminate theory (CLT) was used [3] 
and transverse deformations through the wall thickness were neglected to make the contour 
analysis easier. In advanced aircraft applications, the thin wall assumption of CLT is not 



valid for low aspect ratios which typically use moderately thick-walled sections. For 
composite structure in which strong elastic coupling exists, these transverse stresses and 
strains heavily influence the structural behavior. The effect of transverse shear stress was 
shown to be critical even in the buckling of so-called “ thin ” laminates in a study conducted 
by Chattopadhyay et al [4].. Secondly, in some of the work, the cross sectional geometry 
was assumed to remain rigid during beam deformation and thus in-plane warping was 
neglected [1-2]. However, in-plane warping is important for loaded wing structures with 
short aspect ratio. Studies by Weisshaar et al. [5] indicate that chordwise bending mode, 
which is associated with in-plane warping, is important in the prediction of aeroelastic 
performance. Therefor, in the present research, a composite box beam model is developed 
based on a higher-order laminate theory, that can effectively predict the dynamic response 
under unsteady aerodynamic loads. This model accounts for through-the-thickness 
variations in shear strains and includes both inplane and out-of-plane warping deformation. 

Composite Modeling 

A rectangular composite box beam model with taper and sweep is developed to 
represent the load carrying member of an aircraft wing (Fig. 1). The single-celled 
composite box beam model is based on a higher-order composite laminate theory [6] and 
accounts for the distributions of shear strains through the thickness of each wall. The 
displacement field for each wall section is described by bending, warping and inplane 
stretching. Continuity between the wall displacement fields is imposed at each of the four 
comers of the cross section. The analysis is capable of modeling low aspect ratio wings 
with moderately thick-walled load carrying members while accounting for chordwise 
bending during flutter analysis. The finite element method is used to formulate the 
governing equations of motion. 

For each of the individual plates, the higher-order displacement field is defined in 
local coordinate system as follows (Fig.l). 

u(x,y,t) = u 0 (x,y,t) + zy x (x,y,t) + z 2 £ x (x,y,t) + z 3 £ x (x,y,t) 

v(x,y,t) = v 0 (x,y,t) + z\|/ y (x,y,t) + z 2 ^(x,y,t) + z 3 ^ y (x,y,t) (1) 

w(x,y,t) = w 0 (x,y,t) 

where u 0 , v 0 and w 0 denote the displacements of a point (x,y) on the midplane and \|/ x 
and \|/ y are the rotations of the normal to the midplane about the y and x axes, respectively. 
The higher-order terms £ x , C x , £ y and £ y represent beam warping in each plate. By 
imposing the necessary conditions that the transverse shear stresses must vanish on the 
plate top and bottom surfaces, the following refined displacement field is obtained. 
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u = u 0 + z 


u = u 0 +z 

W = Wq 


4z 2 ,<9w 

y ' _ 3h r( * +y * ) . 
4z 2 .dw 

y * _ 3h r< * +y * ) 


( 2 ) 


where h is the plate thickness. Making the assumption of small displacements and 
rotations, a linear strain-displacement relationship is used. The governing equations of 
motion for an individual plate is derived using the Hamilton’s principle. 

J' | ! S[U,+V,-W,]dt=0 (3) 

where U, V and W denote the kinetic energy, the strain energy and the work done by 
external forces, respectively. Using the constitutive relations for general orthotropic 
material along with the strain-displacement relations, the element stiffness matrix, the mass 
matrix and the forcing vector are derived from Eqn. 3. 

The construction of the box beam from plate elements is shown in Fig. 2. The 
quantities u, v, w are displacements along x, y and z axis, respectively, and 0 X , 0 y and 0 7 

are rotations along theses directions. To make stiffness transformation possible, continuity 
of displacements and rotations are imposed at each of the four comers while the generalized 
forces corresponding to higher order warping terms are set to zero. Through the use of 
coordinate transformation, the reduced stiffness matrix is expressed in the global form. 
Assembly of the element matrices leads to the following governing equation for wing 
dynamic motion. 

Mx + Cx + Kx = q (4) 

where M , C and K denote the global mass, damping and stiffness matrix, respectively. 
The vector x represents structural elastic deformation and the vector q denotes forcing 

vector. 


Results 

The static results of the composite box beam analysis are presented in this section. 
The correlations are made with experimental data which was furnished by studies 
conducted by Chandra et al. [7]. The test beams are single-celled rectangular box beams 
with three kinds of stacking sequences: cross-ply, symmetric and antisymmetric lay-up. 
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The beam dimensions are shown in Fig. 3. Table 1 presents the elastic properties for the 
composite material which is used to fabricate these beams. 

Using the higher-order theory, the static response of the composite box beams is 
computed under unit tip bending and unit tip torque load. Results are shown in Figs. 4-10. 
It is shown that, for cross-ply and symmetric lay-up cases, correlation between higher- 
order solution and test data are excellent. Significant differences exist between 
experimental data and computational results for the (0/30)3 antisymmetric lay-up case. For 
beams with antisymmetric lay-up under unit tip torque, comparisons are also made between 
the deflections predicted by higher-order theory and computations using three dimensional 
finite element formulation. Both solutions give identical results. 

Future Plan 

Following is an outline of the future research. 

(1) Extend the above analysis technique to model two-cell composite box beams. 

(2) Couple the developed composite analysis procedure with aeroelastic technique. 

(3) Perform aeroelastic tailoring of composite airplane wing using a hybrid optimization 
method. 

Briefing describe of the proposed research follows. 

The higher-order box beam theory will be extended to model multi-cell box beam 
configurations in order to deal with more realistic wing structural layouts. This capability 
will be verified by applying the analysis procedure to a two-cell composite box beam 
model. 

Next, the higher-order composite box beam theory will be coupled with unsteady 
aerodynamic computation to perform aeroelastic analysis. Aerodynamic loads will be 
computed using the Doublet Lattice Method (DLM) as implemented in the analysis code 
ASTROS [8]. Both static and dynamic aeroelastic analysis, including predictions of flutter 
and divergence boundaries, will be conducted in the frequency domain. Using the modal 
approach, problems will be solved in the reduced space. This will substantially save 
computational time without compromising the accuracy of the results. Thus the procedure 
will be more suitable for design optimization studies. 

Finally, aeroelastic tailoring will be performed based on the developed higher-order 
composite laminate theory and the aeroelastic analysis. The complex composite tailoring 
problem will be addressed using a hybrid optimization technique. The hybrid technique 
will allow the inclusion of both continuous design variables, such as the wing planform 
geometry, and discrete variables, such as stacking sequence, in the investigation. 
Objective functions and constraints pertaining to improvement in both structural and 
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aeroelastic performance will be included using formal multiobjective function formulation 
technique developed by Chattopadhyay and McCarthy [9]. 
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L/d=56 L/d=29 


L (in) 30 30 

d (in) 0.537 1.025 

c (in) 0.953 2.060 

Ply thickness (in) 0.005 0.005 

Wall thickness (in) 0.03 0.03 


Figure 3 Test beam dimensions 


Table 1. Material properties. 


E l { [ msi) 

20.59 

E T (msi) 

1.42 

G lt ( msi) 

0.87 

G-jji msi) 

0.5 

y it 

0.42 

Ytt 

0.42 


7 






Bending slope (rad) 



Figure 4. Bending slop under unit tip bending load of cross-ply lay-up beam; 
(0/90)3, L/d=29. 
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Figure 5. Twist under unit tip torque of cross-ply lay-up beam; 
(0/90)3, L/d=29. 


Figure 6. Twist at x/L=0.5 under tip torque of anti-symmetric lay-up beams; 
L/d=56. 





Bending slope (rad) 



Figure 7. Bending slope under unit tip bending load of symmetric lay-up beam; 
top & bottom ( 15 ) 6 , sides ( 15 /- 15)3, L/d=56. 
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Twist angle (rad) 



Figure 8. Twist under unit tip bending load of symmetric lay-up beam; 
top & bottom ( 15 ) 6 , sides ( 15 /- 15)3, L/d=56. 
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Bending slope (rad) 



Figure 9. Bending slope under unit tip bending load of symmetric lay-up beam; 
top & bottom (30)6, sides (307-30)3, L/d=56. 
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Twist angle (rad) 



Figure 10. Twist under unit tip bending load of symmetric lay-up beam; 
top & bottom (30)6, sides (30/-30)3, L/d=56. 
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